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Abstract 

DNA methylation (addition of methyl groups to cytosines which normally represses gene 
transcription) and changes in telomere length (TTAGGG repeats on the ends of chromosomes) 
are two molecular modifications that result from stress and could contribute to the long-term 
effects of intrauterine exposure to maternal stress on offspring behavioral outcomes. Here, we 
measured methylation of Brain-derived neurotrophic factor (Bdnf), a gene important in 
development and plasticity, and telomere length in the brains of adult rat male and female 
offspring whose mothers were exposed to unpredictable and variable stressors throughout 
gestation. Males exposed to prenatal stress had greater methylation (Bdnf IV) in the medial 
prefrontal cortex (mPFC) compared to non-stressed controls. Further, prenatally-stressed males 
had shorter telomeres than controls in the mPFC. This study provides the first evidence in a 
rodent model of an association between prenatal stress exposure and subsequent shorter brain 
telomere length. Together findings indicate a long-term impact of prenatal stress on DNA 
methylation and telomere biology with relevance for behavioral and health outcomes, and 
contribute to a growing literature linking stress to intergenerational epigenetic alterations and 


changes in telomere length. 
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1. Introduction 


Perturbations during gestation in the form of malnutrition, infection, or psychosocial 
stress have devastating consequences on the physical and mental health well-being of offspring. 
For example, people exposed to famine during gestation (such as that of the Dutch Famine of 
1944) have significant physical health problems including cardiovascular disease and high blood 
pressure (Painter et al., 2008; Roseboom et al., 2011). Gestational famine exposure also produces 
long-lasting neurobiological and behavioral consequences, including an increased incidence of 
schizophrenia (Hoek et al., 1996; Roseboom et al., 2011). Psychosocial stress experienced by the 
mother during gestation likewise produces vulnerability to psychopathology in offspring, 
including increased risk for schizophrenia and depression (Betts et al., 2015; Malaspina et al., 
2008). Further, work from various animal models provides additional evidence that a suboptimal 
intrauterine environment has long-term physical and behavioral consequences for offspring (for 
reviews see (Bock et al., 2015; Driscoll and Barr, 2016; Lupien et al., 2009)). 

One way that prenatal stress could affect brain function and behavior is through 
epigenetic modifications, including DNA methylation. DNA methylation is characterized by the 
addition of methyl groups to cytosines, typically of cytosine-guanine dinucleotides (CG sites), 
which can dampen transcriptional activity through the blockade of transcription factors or the 
recruitment of transcriptional silencing machinery. Recent studies have also highlighted the role 
of active DNA demethylation in regulating transcription within the CNS, and candidates thought 
to reverse DNA methylation include ten-eleven-translocation (TET) proteins (Kohli and Zhang, 
2013), GADD45b (Ma et al., 2009), and MBD2 (Detich et al., 2002). Though traditionally 
attributed to early cell development and differentiation, there is now a vast literature showing 
that DNA methylation/demethylation remains an active process even occurring in mature 


neurons. Further, studies in humans (Cao-Lei et al., 2014; Kundakovic and Champagne, 2011; 


Nye et al., 2014) and animals (Boersma et al., 2014; Dong et al., 2015; Jensen Pefia et al., 2012; 
Matrisciano et al., 2013; Morgan and Bale, 2011; Mueller and Bale, 2007, 2008) have 
highlighted the ability of prenatal stress to alter DNA methylation at various gene loci, including 
the Brain-derived neurotrophic factor (Bdnf) gene. Bdnf is crucial for development and plasticity 
in the developing and adult brain and is implicated in various psychiatric disorders, including 
depression, schizophrenia, and PTSD (Martinowich et al., 2007). Given the widespread 
functioning of Bdnf across brain regions throughout early development, this locus may play a 


role in translating prenatal stress to adult phenotypic outcomes including psychopathology. 


Recent work has also sought to characterize the detrimental impact of early-life stress by 
examining changes in telomere length —stretches of TTAGGG repeats located on the ends of 
chromosomes (for a comprehensive review, see (Blaze et al., 2015)). Studies continue to show 
that various forms of early-life stress (e.g., perinatal complications, psychological stress, physical 
abuse, etc.) during the prenatal (Entringer et al., 2011; Entringer et al., 2013; Shalev et al., 
2014a) or postnatal (Asok et al., 2013; Kananen et al., 2010; Kiecolt-Glaser et al., 2011; 
O'Donovan et al., 2011; Tyrka et al., 2010) period shorten telomere length in the periphery. 
Given that telomeres are dynamically regulated during the first years of life (Slagboom et al., 
1994; Zeichner et al., 1999), examining neural changes in telomere length following prenatal 
stress is another important avenue for understanding how stress in early-life may relate to future 


mental well-being. 


In daily life, humans tend to be exposed to unpredictable challenges that vary in 
magnitude and duration. Studies of prenatal stress in rodents often employ repeated presentations 
of a single, non-varying stressor, most frequently, immobilization stress (de Souza et al., 2013; 


Jensen Pefia et al., 2012; Matrisciano et al., 2013). However rodent studies provide evidence for 


habituation to predictable, non-varying stressors (Armario, 2006; Daviu et al., 2014), and studies 
employing variable prenatal stress paradigms have reported etiological changes that may be more 
closely linked to psychiatric disorders (Bock et al., 2015; Mueller and Bale, 2008; Sickmann et 
al., 2015; Wilson et al., 2012). Collectively, these findings suggest that variable, unpredictable 


stress may model the human condition with greater fidelity than invariant stress paradigms. 


The current study was designed to characterize changes in Bdnf DNA methylation and 
telomere length across multiple discrete brain regions of the adult rat following unpredictable 
variable prenatal stress (UVPS). We chose to focus on the medial prefrontal cortex, ventral 
hippocampus, dorsal hippocampus, and central/basolateral amygdala, as our prior work showed 
changes in Bdnf exon IV DNA methylation and telomere length in these regions in adult animals 
subjected to maternal maltreatment (Asok et al., 2014; Blaze et al., 2013; Roth et al., 2009; Roth 
et al., 2014). Changes in Bdnf 1V DNA methylation have likewise been found in these same 


brain regions in other models of prenatal stress (Boersma et al., 2014; Dong et al., 2015). 


2. Materials and Methods 


2.1 Subjects 


Female Sprague-Dawley (SD) rats and their offspring were used (n=6 control animals per 
sex, 7 UVPS animals per sex derived from 13 dams, with only one male and one female 
offspring from each control and UVPS litter). To avoid stress associated with shipping during 
pregnancy, dams were time-mated in our breeding colony (Taconic Farms, Germantown, NY). 
Pregnancy was determined by daily vaginal lavage with the presence of spermatozoa indicating 


Gestational day zero (GO). Throughout pregnancy, dams were housed in groups of three in 


standard maternity cages (47cm x 26cm x 21cm) lined with corncob bedding and maintained 
under standard colony conditions (12:12 light/dark cycle [0600:1800]; 21 +/- 1° C at 30 - 50% 
humidity). Rat chow (Purina #5102) and water were available ad libitum. On gestational day (G) 
21, dams were singly housed. Animal experimentation was conducted in full accordance with the 
NRC Guide for the Care and Use of Laboratory Animals (National Academy of Sciences, 
copyright 2011). All procedures were approved by the Wake Forest School of Medicine Animal 


Care and Use Committee (ACUC). 


2.2 Unpredictable, Variable Prenatal Stress (UVPS) Paradigm 


On GO, pregnant dams were randomly assigned to either control or UVPS treatment conditions 
(Body weight: CNTL, 288 + 26g; UVPS, 292 + 23g). On G1, UVPS dams were exposed daily to 
three stressors, one during each of three discrete time periods varying in duration (Time periods, 
Early: 0600-1200, Mid: 1200-1800, Late: 1800-600; Duration, 15, 30 or 60 minutes). Stressors 
included tube restraint (PVC tube, length x diameter: 21.6 x 6.4 cm), high-frequency strobe light 
(60Hz, 3.5W, 85mA), and white-noise induced stress (80db). Stimulus presentation was block 
randomized across pregnancy and spanned G1 — G21. Control dams were handled briefly each 
day to match UVPS handling. To eliminate potential confounds of altered postnatal maternal 
care by UVPS dams contributions on experimental endpoints, offspring in both conditions were 


fostered at birth to non-manipulated, newly parturient dams. 


2.3 Tissue Collection and Processing 


On postnatal day (PN) 90, rats were euthanized and flash frozen. Brain were extracted 
and later sectioned into 300 um slices on a cryostat. DNA and RNA were simultaneously 
extracted from homogenized tissue punches of medial prefrontal cortex (mPFC; consisting of 
infralimbic and prelimbic cortices), hippocampus (dorsal vs. ventral) and central/basolateral 
amygdala (Qiagen Inc., Valencia, CA). Quantity and quality of nucleic acids were determined 


using spectrophotometry (Nanodrop 2000) and then frozen at -80°C for later processing. 


2.4 Bdnf DNA Methylation and Gene Expression Assays 


For methylation assays, DNA was bisulfite-converted (Epitect Bisulfite Kit, Qiagen, Inc., 
Valencia, CA) and direct bisulfite-sequencing PCR (BSP) was used to measure methylation at 11 
CG sites in Bdnf exon IV as previously reported (Roth et al., 2009; Roth et al., 2014). PCR 
products were purified (Diffininty Genomics, RapidTip) and sequenced using reverse primers at 
the Delaware Biotechnology Institute. The electropherogram was read on Chromas software, 
where the percent methylation of each CG site was determined by the ratio between peak values 
of G and A (G/[G +A]). We confirmed the accuracy of this technique by assessing methylated 
standards (Rat Pre-mixed Calibration Standards, Epigendx) ranging from 0-100% methylation 


and performing a linear regression on methylation values (F (1,40)=228.9, p<0.001, r7=0.8512). 


For gene expression assays, RNA was reverse transcribed using a cDNA synthesis kit 
(Qiagen, Inc., Valencia, CA) and amplified via real-time PCR using Taqman probes 
(ThermoFisher Scientific, Grand Island, NY) to target Bdnf exon IV-containing transcripts. 


Tubulin was used as a reference gene and all reactions were run in triplicate. The comparative Ct 


method was used to obtain the relative Bdnf 1V mRNA fold change in experimental (UVPS) 


relative to control animals (Livak and Schmittgen, 2001). 


2.5 Telomere Length and Telomerase assays 


DNA was assessed for purity, and diluted to 10ng/uL via Nanodrop spectrophotometry 
following procedures used previously(Asok et al., 2014). Quantitative real-time PCR was 
conducted with the following primers targeting telomeres (T; Forwardret:5’- 
CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT-3’ and ReversersgL:5’- 
GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT-3’) and a single copy control 
gene (S), the angiotensin receptor gene NCBI M86912.1 (Forwardarti: 5’- 
ACGTGTCTCAGCATCGACCGCTACC -3’ and Reverseari: 5’- 
AGAATGATAAGGAAAGGGAACAAGAAGCCC -3’) (O'Callaghan and Fenech, 2011). Each 
PCR well contained a final concentration of 1x Power Sybr Green Master Mix (ThermoFisher 
Scientific, Grand Island, NY), 100 nM forward primer, 100 nM reverse primer, and 20ng of 
sample DNA. The experimenter was blind to a sample’s group identification when loading the 
PCR plate. For each sample, the telomere and AT1 qPCR assays were carried out in triplicate in 
the same well position on different 96-well plates. All male and female samples for a particular 
brain region were run on the same 96-well PCR plate. Given that plates did not contain a 
standard, we did not statistically contrast brain regions across PCR plates. All samples were 
assayed in triplicate and any replicate that deviated + | cycle threshold (Ct) beyond the triplicate 
average (~ 4.3% of all replicates) was excluded from the final Ct average. Eight samples failed to 


amplify in the amygdala (2 females and 2 males in the control condition and | female and 3 


males in the UVPS condition) and were thus excluded from analyses. Relative telomere length 
was calculated as a ratio of telomeres (T) to single copy gene (S) by the formula T/S = (24* 


tel)(24C‘AT!). The T/S ratio was then transformed into a proportionate score against the 


unstressed control condition for each brain region. 


Similar to Bdnf IV gene expression assays, CDNA was amplified via real-time PCR using 
Taqman probes (ThermoFisher Scientific, Grand Island, NY) that targeted telomerase exons 7-9- 
containing transcripts. Tubulin was also used as a reference gene and all reactions were run in 
duplicate. The comparative Ct method was used to obtain the relative mRNA fold change in 


experimental (UVPS) relative to control animals (Livak and Schmittgen, 2001). 


2.6 Statistical Analyses 


For each brain region, BSP data were analyzed with two-way ANOVAs comparing 
average methylation (across all 11 sites) between treatment group and sex (2 x 2), as well as two- 
way ANOVAs on each of the 11 CG sites [treatment group (2) x sex (2)]. Telomere length 
between conditions and sexes was likewise analyzed with two-way ANOVAs for each brain 
region. Bonferroni multiple comparison tests were used when appropriate for post hoc analyses. 
Differences in Bdnf or telomerase mRNA levels were analyzed by one-sample t-tests (for 
comparison of UVPS subjects to controls). Differences were considered to be statistically 


significant for p<0.05, and non-significant trends at p<0.1 are also reported. 


3. Results 


3.1. Bdnf IV DNA methylation 


To assess the effects of UVPS on adult DNA methylation patterns, we measured 
methylation of Bdnf IV in adult (PN90) male and female rats. The primer set we used targeted 
DNA that encompasses a transcription factor binding and transcription start site. We targeted this 
region as there are ample data indicating that this region of the Bdnf gene seems especially 
sensitive to a range of environmental (physical or psychosocial) factors (Kundakovic et al., 2013; 


Lubin et al., 2008; Roth et al., 2009; Roth et al., 2014; Roth et al., 2011). 


Medial prefrontal cortex (nPFC). UVPS males had higher levels of methylation in the 
mPFC (Figure 1A). A two-way ANOVA comparing average methylation values for both 
treatment (UVPS vs. control) and sex (males vs. females) revealed a significant interaction 
(F(1,19)=18.19, p=0.0004). UVPS males had higher overall methylation than control males 
(t=3.512, p=0.0140) and UVPS females (t=3.374, p=0.0191). UVPS males also had higher 
methylation at several individual CG sites. A two-way ANOVA comparing methylation values at 
CG site 1 for both treatment (UVPS vs. control) and sex (males vs. females) revealed a 
significant interaction (F(1,19)=16.66, p=0.0006). UVPS males had higher methylation than 
control males (t=3.165, p=0.0306) and UVPS females (t=3.552, p=0.0128). A two-way ANOVA 
comparing methylation values at CG site 2 for both treatment and sex revealed a significant 
interaction (F(1,19)=12.31, p=0.0023). Control females had marginally higher methylation than 
UVPS females (t=2.854, p=0.0609) and control males (t=2.794, p=0.0694). A two-way ANOVA 
comparing methylation values at CG site 3 for both treatment and sex revealed a significant 


interaction (F(1,19)=16.75, p=0.0006). UVPS males had higher methylation than control males 
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(t=2.921, p=0.0525) and UVPS females (t=3.366, p=0.0195). A two-way ANOVA comparing 
methylation values at CG site 4 for both treatment and sex revealed a significant main effect of 
treatment (F(1,17)=13.75 p=0.0017), sex (F(1,17)=6.453, p=0.0211), and an interaction 
(F(1,17)=25.35, p=0.0001). UVPS males had higher methylation than control males (t=6.052, 
p<0.0001), control females (t=4.325, p=0.0028), and UVPS females (t=5.477, p=0.0002). 
Finally, a two-way ANOVA comparing methylation values at CG site 7 for both treatment and 
sex revealed a significant interaction (F(1,19)=14.38, p=0.0012). UVPS males had higher 


methylation than control males (t=3.452, p=0.0160) and UVPS females (t=3.145, p=0.0320). 


Ventral hippocampus. As depicted in Figure 1B UVPS females had somewhat higher 
levels of methylation for several CG sites, though most statistical comparisons failed to reach 
significance. A two-way ANOVA comparing average methylation values for both treatment 
(UVPS vs. control) and sex (males vs. females) revealed a marginally significant interaction with 
females showing higher overall methylation (F(1,22)=3.917, p=0.0605). A two-way ANOVA 
comparing methylation values at individual sites for both treatment (UVPS vs. control) and sex 
(males vs. females) revealed a significant interaction at CG site 1 (F(1,22)=5.672, p=0.0263) and 
CG site 2 (F(1,21)=4.935, p=0.0374), though all post hoc analyses failed to reach statistical 
significance (for site 1 however the comparison between control and UVPS females reached 


marginal significance (t=2.607, p=0.0966)). 


Dorsal hippocampus and amygdala. There were no main effects of treatment or sex, nor 
any interaction effects for average methylation levels or methylation at individual CG sites 


within the dorsal hippocampus or central/basolateral amygdala. 
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3.2 Bdnf IV gene expression 


As one way to envision the functional relevance of our changes in DNA methylation 
induced by UVPS within the mPFC, we measured steady-state gene expression of Bdnf exon IV- 
containing transcripts. No significant differences in Bdnf IV mRNA levels were detected 
between UVPS males and controls (t(6)=1.900, p=0.1061) nor UVPS females and controls 
(t(6)=1.693, p=0.1415). These data however do not necessarily mean our observed methylation 
changes are not functionally relevant in terms of Bdnf gene expression. DNA methylation 
changes do not always result in steady-state changes in gene expression, but can instead prime 
transcriptional responses to subsequent stimuli and neural activation (Baker-Andresen et al., 


2013). 


3.3. Telomere length 


Medial prefrontal cortex (mPFC) telomere length. UVPS had a significant effect on 
telomere length in the mPFC (Figure 2). A two-way ANOVA revealed a main effect of treatment 
(UVPS vs. control) on telomere length in the mPFC (F(1,22)=7.847, p=0.0104). UVPS males 
exhibited shorter telomeres than control males (t(11)=2.230, p=0.048). In females, UVPS 
females exhibited marginally shorter telomeres relative to controls (t(11)=1.902, p = 0.084). 

Ventral hippocampus, dorsal hippocampus, and amygdala telomere length. In the 
hippocampus and central/basolateral amygdala there were no main effects of treatment or sex, 


nor any interaction effects on telomere length. 
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3.4 Medial prefrontal cortex (mPFC) telomerase mRNA 

Given the differences observed in telomere length within the mPFC, we also wanted to 
selectively examine telomerase in this brain region. No significant differences in telomerase 
were detected between UVPS males vs. controls (t(6)=1.807, p=0.120), but there was a marginal 
decrease in telomerase mRNA for UVPS females compared to control females (t(6)=2.265, 


p=0.064). 


4. Discussion 


Here we used a rodent model of unpredictable, variable prenatal stress to test the 
hypothesis that a stressful environment in utero can render epigenetic and telomeric changes in 
the brain of offspring. We observed that Bdnf 1V DNA methylation patterns in offspring at age 
PN90 varied depending on prenatal condition, sex, and brain region. Specifically, males from 
prenatally stressed mothers had significantly higher levels of methylation in the mPFC compared 
to controls. This effect was not present in females from stressed mothers, who instead had lower 
levels of methylation than prenatally stressed males within the mPFC and marginally higher 
levels of methylation than controls within the ventral hippocampus. No differences in Bdnf IV 
DNA methylation were found in the amygdala or dorsal hippocampus in the stressed group 
compared to controls for either sex. Similarly brain-region specific, stress-induced effects on 
telomere length were detected — male, and to a lesser extent female, UVPS offspring exhibited 


shorter telomeres in the mPFC with no differences in the amygdala or dorsal hippocampus. 


While decades of research have shown the detrimental effects of prenatal stress on brain 
function and behavior, there remains much to elucidate regarding the molecular mechanisms 


underlying these effects. A large drawback of this field as a whole is the variance in prenatal 
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stress paradigms used, which makes direct comparisons between studies difficult (Bock et al., 
2015). The timing of gestational stress, predictability of stressor, and the stressors used vary 
greatly between studies, but are extremely important factors when considering the 


neurobiological changes induced by prenatal stress. 


Studying PN21 and PN80 male offspring of pregnant rats subjected to a different variable 
stress paradigm (beginning on G14, dams were subjected to a mixture of restraint stress, forced 
swim, cold exposure, social stress, and light exposure overnight) than the one we employed here, 
another lab found that prenatally stressed males at both ages had higher Bdnf 1V DNA 
methylation in the amygdala and hippocampus (only a site-specific change at PN80) (Boersma et 
al., 2014). One CG site showed higher methylation within the PFC of prenatally stressed males 
at PN21, however no methylation differences were detected in the PFC of PN80 offspring 
(Boersma et al., 2014). Studying PN75 male offspring of pregnant mice subjected to restraint 
stress under bright light conditions (45 minutes three times per day from G7 through delivery), 
another lab observed significantly higher levels of methylation (as well as hydroxymethylation) 
for DNA associated with Bdnf IV in the frontal cortex (Dong et al., 2015). Finally, studying 
PN90 male and female offspring of pregnant mice exposed to predator stress (from G11-18 dams 
were daily exposed to three different predator urines), another lab found lower methylation of 
Bdnf IV DNA (at only one CG site) in the hippocampus of female offspring (St-Cyr and 


McGowan, 2015). 


Differences between these studies and ours may be attributable to differences in the type 
and timing of prenatal stress. Our study is unique in that we employed a block randomized 
program of three stressors (light, sound and restraint) presented early- and mid-day during the 


lights on phase of the cycle and once during the dark phase for varying durations ranging from 
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15 to 60 minutes for the entire span of gestation. In addition, examining methylation within the 
entire prefrontal cortex, amygdala, and hippocampus instead of the specific divisions within 
could produce different results. Nonetheless, data altogether are consistent with the notion that 
stress experienced by the dam is able to produce long-term epigenetic changes in offspring at a 
gene well-linked to neuroplasticity and cognition (e.g.(Bekinschtein et al., 2008; Dincheva et al., 
2016; Martinowich et al., 2007; McEwen et al., 2015)) and mental health (e.g.(Fuchikami et al., 
2011; Keller et al., 2010; Perroud et al., 2013; Smith et al., 2011)). Our study and that from the 
McGowan lab shed light on the sex-specific epigenetic changes inducible by intrauterine 


exposure to maternal stress. 


Here we also found that prenatally stressed males exhibited shorter telomeres in the 
medial prefrontal cortex compared to controls. A number of studies have demonstrated a strong 
relationship between early-life stress, reduced peripheral telomere length, and future mental- 
health (Entringer et al., 2011; Marchetto et al., 2016; Shalev et al., 2014b). The effects of stress 
on telomere length in the brain is less clear (Thomas et al., 2008). Telomeres shorten with age in 
the cortex of rodents, but this shortening is thought to reflect changes occurring within glial cells 
given their increased mitotic potential relative to neurons (Flanary and Streit, 2003; Linkus et al., 
2016; Szebeni et al., 2014). We did not distinguish between neurons and glia in measuring 
telomere length or methylation in the present study, although recent studies from our lab suggest 
that developmental (postnatal) stress induces changes in mPFC Bdnf IV DNA methylation 


primarily in neurons (Blaze and Roth, 2017). 


Rats selectively bred for enhanced depressive-like symptoms (i.e., the Flinders-sensitive 
line) show reduced telomere length, telomerase, and Bdnf expression in the hippocampus (Wei et 


al., 2015) - supporting a role for telomere shortening in the brain. In contrast, other studies have 
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found that postnatal stress (e.g., maternal separation or caregiver maltreatment) is associated 
with longer telomeres in the ventral hippocampus and medial prefrontal cortex (Asok et al., 
2014; Botha et al., 2012). The findings here complement our previous work with regard to stress 
effects on a common brain region (i.e., the medial prefrontal cortex), but differ in the direction of 
telomere change. Similar to methylation, it is likely that the type (maltreatment vs. maternal 
separation vs. variable stress) and timing (pre- vs. postnatal) of stress can differentially affect 
brain telomere length, and future studies are certainly warranted to better understand the 


parameters affecting brain telomere length. 


In summary our findings highlight the capacity of in utero stress to leave its epigenetic 
and telomeric mark within specific regions of the brain, providing important new insights into 
how gestational stress may have far reaching behavioral consequences for progeny. Indeed, 
behavioral analysis of the prenatally stressed and non-stressed subjects used in this report 
revealed robust sex-specific (male-biased) effects on anxiety responses (Ronca, unpublished). 
These findings are consistent with reports that prenatal stress leads to heightened stress 
sensitivity, especially among males (Bronson and Bale, 2014), including studies employing 
variable prenatal stress paradigms (Koenig et al., 2005; Wilson et al., 2013). In accord with 
previous research linking prenatal stress with changes in body weight regulation and metabolic 
disorders (Boersma et al., 2014; Wilson et al., 2013), adult body weights of prenatally stressed 
male (but not female) offspring of the present study were significantly elevated relative to non- 
stressed offspring (Ronca, unpublished). Our findings indicate that, in addition to brain 
epigenetic and telomeric changes, unpredictable, variable gestational stress exerts clear sex- 


specific programming effects on the physical and behavioral health of the offspring. 
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The study reported herein provides the first evidence in a rodent model for an enduring 
influence of prenatal stress on both DNA methylation and telomere length in adulthood. Going 
forward, the challenge will be to identify the specific role(s) of these molecular brain changes in 
relation to the diverse effects of early life stressors. Recent progress characterizing physiological 
and structural changes in the brain associated with prenatal stress, and their sex-specificity, may 
provide important new insights. Neuronal changes are highly region-specific with the most 
dramatic changes occurring in limbic and prefrontal cortical areas, those involved in cognition 
and emotional behavior. Altered neurogenesis, neuronal arborization, neuronal density, dendritic 
architecture, synaptic connectivity, and myelination have all been reported in animal models of 
prenatal stress (Bock et al., 2015; Schuurmans and Kurrasch, 2013; Soztutar et al., 2016). 
Telomere length may be related to changes in a variety of cells throughout the brain, including 
glia (Dabouras et al., 2004; Leventopoulos et al., 2007). Relating these cellular alterations to 
epigenetic and telomeric changes in common brain regions will likely lead to deeper 


understanding of the underlying mechanisms. 


While extensive progress has been made in the understanding of the role of early life 
adversity in the development of psychiatric and metabolic disorders, more detailed research on 
sex-specific epigenetic brain structural and behavioral changes is warranted. In particular, 
consistency across gestational stress paradigms, including the predictability and timing of 
stressor, the use of common behavioral and endpoint measures, and comparisons between male 
and females are vital for uncovering consistent and meaningful patterns of neurobiological 
changes induced by prenatal stress. Regardless, our present findings have strong relevance for 
behavioral and health outcomes, and contribute to a growing literature linking stress to early 


developmental programming, intergenerational epigenetic transmission of environmental 
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exposures, and telomere biology. Illuminating molecular mechanisms underlying 
transgenerational programming of stress responses and mental and physical health outcomes are 
important first steps toward establishing therapeutic interventions to halt the generational 


persistence of pathological outcomes associated with early life adversity. 
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Figure Captions 


Figure 1. Bdnf IV DNA methylation in the A) medial prefrontal cortex and B) ventral 
hippocampus. UVPS male offspring had significantly higher methylation compared to control 
males and UVPS females across most CG sites within the medial prefrontal cortex. *p<0.05, 


**H<0.01, ***p<0.001; n=6-7/sex/condition; error bars represent SEM. 


Figure 2. Telomere length in the adult medial prefrontal cortex. Males who experienced UVPS 
had significantly shorter telomeres compared to control males. UVPS females also exhibited 
marginally shorter telomeres relative to controls. *p<0.05; n=6-7/sex/condition; error bars 


represent SEM. 
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